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Abstract 

A thermodynamic foundation using the concept of internal state variables is presented for the 
kinematic description of a viscoplastic material. Three different evolution equations for the 
back stress are considered. The first is that of classical, nonlinear, kinematic hardening. The 
other two include a contribution that is linear in stress rate. Choosing an appropriate change in 
variables can remove this stress rate dependence. As a result, one of these two models is shown 
to be equivalent to the classical, nonlinear, kinematic hardening model; while the other is a new 
model — one which seems to have favorable characteristics for representing ratchetting behavior. 

All three models are thermodynamically admissible. 

Let us consider a complementary (Gibbs) free energy for the thermodynamic potential function, 

i.e. 

= y[T,(j uv ,f3 uv \ , Pu = 0 , ( 1 ) 

where square brackets [•] are used to denote ‘function of’. The temperature T and stress 
are external variables whose variations can, in principle, be controlled by an observer. However, 
the evolution of the internal state variable flij — associated with kinematic hardening — cannot be 
controlled by an observer; it is a material response. This variable must therefore evolve as a function 
of state. Furthermore, if this kinematic variable is to remain invariant under transformation to 
another kinematic variable, then its evolution will, in general, also depend upon the rates of change 
in the external variables [1]. For simplicity, let us consider the following equation for the evolution 
of internal state, i.e . 

Pi 3 = Xu P\ ^ , M + (N/2H) Sij , (2) 

where Sij = an - 1 k cr kk^ij is the deviatoric stress with S tJ representing the Kronecker delta, and 
where N and H are material constants. In general, a term that is linear in T may also be introduced 
into (2) [2], but it is ignored in this paper for the purpose of simplification. Notice that /3;j, and 



therefore x«j, are deviatoric by definition. The function Xij describes the irreversible contribution 
to the evolution of 

Given (1) and (2), one can derive [1, 2] relationships that govern the entropy, 


s = 


-dV 
dT ’ 


(3) 


the strain, 


Sij 


and the intrinsic dissipation, 
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(5) 


where e v ^ is the plastic strain, which is deviatoric. It evolves according to a separate evolution 
equation (13). 

If we assume that the complementary free energy (1) is given by 
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thermoelastic contribution 

/ * N 

4^ &ij $ij ~ ~ a AT an — So AT — C |t — To ^1 + In J ^ j + 

Hflijflij + N Sij Sij — Pi j Sij + /3uu]^ » 

^ v / 

viscoplastic contribution 


( 6 ) 


then we obtain — from (3) and (4) — the constitutive equations for an isotropic Hookean material ; 
they are: 

AS = aa ii + (C/T)AT J (7) 

an = 3k (eu - a AT 6a) , (8) 


and 

Si-j = 2/i i^Eij - , 


(9) 


provided that 


OA -N OA 

OSij ~ 2 H dpij ’ 


( 10 ) 


Here /x and k are the elastic shear and bulk moduli, a is the coefficient of thermal expansion, C is 
the specific heat, 5o and To are the initial values of entropy and temperature with AS = S — So 
and AT = T — To, and Eij = eij - £kk6ij is the deviatoric strain. The first three terms in 
the viscoplastic contribution to $ are introduced to remove unwanted cross products that would 
otherwise appear in (9) because of (4) [2]. One also obtains 
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dPij 


= 2 H Pij- 




which defines the thermodynamic force conjugate to /?,j, As a consequence, (5) becomes 


(ii) 


- 2 # Xij 


M\ . 

~Wj) Xi ^ 


( 12 ) 
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which describes the intrinsic dissipation properties of our material. The function A, which is 
constrained by (10), is introduced into $ to affect the intrinsic dissipation (12); its form is model 
dependent. 

For the evolution of plastic strain, we shall consider that 


e v - = V2 II k v II 

*’ /2N "iis-BH 


(13) 


where the back stress, 


Bij = 2 HPh = 2 H X ij + NSij , 


(14) 


accounts for kinematic behavior, with H > 0 being its modulus and 0 < N < 1. The norms used are 
defined by || 1 1| = yj I t] I t] where 7,j is any deviatoric stress-like quantity, and by || J || = yJ2J XJ J i} 
where J,y is any deviatoric strain-like quantity. These von Mises norms are scaled for shear. 


Model I 


The classical, nonlinear, kinematic hardening model [3, 4] has an evolution of internal state de- 
scribed by 

H v ■ 

(15) 


. P Hxa 

Xij = 4 J 


where L is the limiting state for the back stress B{j. There is no stress rate term in the evolution 
of Pij for this model, i.e. N = 0 in (2), and consequently 


B„ = 21 1 (ij - fj- Ill'll) 


(16) 


describes its evolution for the back stress. Because N = 0, the last three terms in the complementary 
free energy (6) do not contribute to the dissipation in this model; its intrinsic dissipation (12) is 
therefore given by 


\S-B\\ + 




( 17 ) 


and it is always satisfied. Hence, Model I is thermodynamically admissible. 


2 Model II 

This model uses the same description for the irreversible evolution of internal state that Model I 
uses, viz . (15). In addition, it assumes that 0 < N < 1, which introduces a reversible attribute to 
the evolution of internal state. As a result, 

B„ = 2 H (if, - B ‘> ~* S ” || i” ||) + NSij ( 18) 

describes this model’s evolution for the back stress. Model II reduces to Model I when N — 0. 

As Lubliner [1] discussed in his paper, one can always (in principle) transform from an internal 
state variable whose evolution contains terms that are linear in the external variable rates, to 
another internal state variable where there are no external variable rates present in the evolution 
equation. This is accomplished in our case by considering the linear transformation 

^ x ' 3 = 2H ' X,} , (19) 
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where the variable Xij is a back stress, but different from 7?,j, with H' = 77/(1 — A) as its associated 
hardening modulus. This transformation enables the flow law (13) to be rewritten in an equivalent 
form as 

* = (20) 
Likewise, it allows the evolution equation for back stress (18) to be rewritten in the equivalent form 


Xij = 2 H' (if. - gf He'll) , 


( 21 ) 


where V = L/( 1 — TV) is the limiting state for the back stress X tJ . Upon comparing (20) and 
(21) with (13) and (16), one observes that they are identical in mathematical structure, but with 
different values for their constants. These differences lead to differences in the intrinsic dissipation 
properties of the material [2]. 

From the perspective of material science [5], the physically correct, internal, state variable has an 
evolution equation with no Sij or T dependence. This coincides with the experimental observation 
that an instantaneous change in either stress or temperature does not produce an instantaneous 
change in a material’s internal structure, t.e. its dislocation structure. Consequently, the back 
stress Xij is the physically correct back stress for Model II, and it is referred to as the physical 
back stress. (The back stress is the physical back stress of Model I.) 

A description of Model IPs dissipation response requires knowledge of its material function A, 
which is found in the expression for the complementary free energy (6). This function is evaluated 
by combining (12), (14), (15), (19), and (20), and then using A to cancel out those terms in 
the dissipation inequality that can become negative valued. This process leads to the differential 


equation 


d\ 

dSij 


N 

1 - N 



( 22 ) 


which is the simplest of several possible solutions. This equation can be integrated, in conjunction 
with the constraint given in (10), to produce 


A = 


N 

1 - N 


N 

4 H 


Sij Sij 



(23) 


which when substituted into (6) defines for Model II; it is the basic constitutive equation of this 
particular model. It follows then that the intrinsic dissipation (12) for Model II is given by the 
inequality 

f||S-X|| + M^ P1>0, (24) 

and it is always satisfied. Hence, Model II is thermodynamically admissible. Notice the similarity 
between (17) and (24). 


3 Model III 

This model considers the irreversible evolution of internal state to be described by 


jr 


le p l 


(25) 
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which differs from (15) by the exchange of Xij with fin in the dynamic recovery term. Like Model II, 
this model takes 0 < N < 1, and therefore 


= 2 H (<?. - I* || i'|l) + NSn 


(26) 


describes the evolution for the back stress Bij. Equation 26 was first proposed by Ramaswamy 
et al. [6]. Notice that the only difference between (16) and (26) for Models I and III is the presence 
of the term N Sij found in (26). Model III reduces to Model I when N = 0. 

Using the same linear transformation that was used in Model II, i.e. (19), one determines that 
Model III has the same transformed flow law (20) as Model II, but it has a different evolution 
equation for the physical back stress Xij, viz . 


- (*p _ V-WXij + NSij 

1 - TV \ ij 2 L 11 


(27) 


To the best of our knowledge, this ia a new expression for the evolution of the physical back stress. 
Here N proportions the dynamic recovery between the physical back stress Xij and the applied 
stress Sij. 

We set out to derive the dissipation response of Model III as we did for Model II. We begin by 
considering a decomposition 

A = A t + A 2 , (28) 

where Ai is taken to be given by (23), i.e. Ai is the A of Model II. The remaining function is 
evaluated by combining (12), (14), (19), (20), and (25), and then using A 2 to cancel out those 
remaining terms in the dissipation inequality that can become negative valued. Because of the 
constraint equation (10), one obtains two partial differential equations, viz. 


<9A 2 

dSij 

<9A 2 

mj 


N 

21(1 - N) 
-H 

L(l-N) 


• SktPkt 

/ S l} - 211 (3 ij 

BPaY 1 

(29) 

Ml S-2H(3\\ 

L ) ’ 

SktPkt ( 

' Sij - 2H(3ij 

nPuY 1 

(30) 

,||S - 2H(3\\ 

L ) ’ 


which when integrated will lead to the complementary free energy $ for Model III, i.e . its fun- 
damental constitutive equation. We have not integrated these equations. For our purpose, it is 
sufficient to know only that A exists. The intrinsic dissipation (12) for Model III is therefore given 
by the inequality 


I5-XII + 


L{ 1 


BJIM 

-N)J 


ll* p ll>o, 


(31) 


and it is always satisfied. Hence, Model III is also thermodynamically admissible. 

A preliminary study of these three models [7] indicates that Model III may be able to predict 
realistic ratchetting behavior; whereas, Model I (and therefore Model II) is known to overpre- 
dict ratchetting behavior [8]. Continued research is required to better understand the predictive 
capabilities of Model III. 
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